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Abstract 

We treat the gravitational effects of quantum stress tensor fluctuations. An operational 
approach is adopted in which these fluctuations produce fluctuations in the focusing of a 
bundle of geodesies. This can be calculated explicitly using the Raychaudhuri equation as 



o 



a Langevin equation. The physical manifestation of these fluctuations are angular blurring 
and luminosity fluctuations of the images of distant sources. We give explicit results for 



the case of a scalar field on a flat background in a thermal state. 
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I. INTRODUCTION 



Although quantum states in field theory may or may not be eigenstates of the 
Hamiltonian, and hence have fixed energy, they are never eigenstates of the energy 
density or other stress tensor components. This means that all quantum field theories 
exhibit fluctuations of the stress tensor. These fluctuations can manifest themselves 
either through fluctuating forces on material bodies, or through fluctuations of the 
gravitational field. In the former category are fluctuations of Casimir forces [l|, H IE HJ 
and fluctuations of the radiation pressure ||. In this paper, we will be concerned 
with the latter phenomenon, where stress tensor fluctuations drive fluctuations of the 
spacetime geometry. These are "passive" fluctuations of the spacetime, in contrast to 
the "active" fluctuations which arise from the quantization of gravity itself. The role 
of stress tensor fluctuations in gravity theory has been discussed by several authors jH, 
0, H II El El El El EI EE EE E3- Among the issues of interest is the extent 
to which these fluctuations place limits on the validity of the semiclassical theory 
of gravity, which is the approximation in which the gravitational field is assumed 
to be classical and determined by the renormalized expectation value of the stress 
tensor [E H EE|. Other authors have been concerned with the possible role of stress 
tensor fluctuations in the early universe or upon black hole evaporation 11 , 0, 

E1E1E3- 

One of the basic problems which must be addressed in treatments of quantum 
stress tensor fluctuations is that the stress tensor correlation function is singular, even 
if the stress tensor operator has been renormalized to have finite expectation values. 
For this reason, many of the earlier treatments @ 0, 0, ©, El studied only a portion 
of the correlation function which is free of singularities, the "fully normal ordered 
part". However, more recently it has become clear that the remaining parts (the 
"vacuum" and "cross" terms to be defined in Sect. [TTJ) can also have physical effects, 
and that extracting these effects requires a careful procedure for defining integrals 
with singular integrands H 0]. For example, the radiation pressure fluctuations, 
arising when light in a coherent state impinges upon a mirror, come entirely from the 
cross term H- 

In this paper, we will present an approach to studying the gravitational effects of 
stress tensor fluctuations which is based upon Raychaudhuri's equation. This equa- 
tion relates the Ricci tensor to the rate of change of the expansion along a congruence 
of geodesies, and plays a central role in the proofs of the classical singularity theorems. 
We will treat it as a Langevin equation which describes how stress tensor fluctua- 
tions lead to fluctuations in the expansion. This in turn leads to a lack of focusing 
and fluctuations in the luminosity of the images of distant objects seen through a 
fluctuating spacetime geometry. Here the quantum matter field will be a massless, 
minimally coupled scalar field on an average background of Minkowski spacetime. 
The outline of the paper is as follows: Section |H] will discuss the stress tensor cor- 
relation function and its decomposition into individual terms. Section IIHI will deal 
with the Ricci tensor correlation function and the fluctuations of the expansion. The 
treatment of the singular pieces in the correlation function will also be discussed in 
this section. In Sec. IIV[ we will explicitly calculate the expansion fluctuations for a 
thermal state. The results of the paper will be summarized and discussed in Sec. fVl 
Except as otherwise noted, we use units in which h = c = G = k=l, where G and 
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k are Newton's and Boltzmann's constants, respectively. 



II. QUANTUM STRESS TENSOR FLUCTUATIONS 

In this paper, we will deal with the fluctuations of the stress tensor for quantum 
fields on an average background of Minkowski spacetime. The formal expectation 
values of T^ v are divergent, so we renormalize by subtracting the expectation value 
in the Minkowski vacuum state: 

'■Tjiv '■— Tjiv — (Tfj,v)o , (1) 

where : T^ v : is the normal ordered stress tensor operator. The fluctuations of : T^ v : 
are described by the stress tensor correlation function 

S^ a p(x,x') = {:T lxv {x)::T aP {x'):) - (:T^(x) :){:T a(5 {x') :) . (2) 

Here the expectation value can be taken in an arbitrary quantum state. It is of 
interest to note that this correlation function is unchanged when : T^ v : is shifted by 
a constant, so we could formally write S^ va p{x,x') without normal ordering: 

S^iuafB (x, X ) = (T„„(x)T a p(x')) - (T, u (x)){T aP (x')) . (3) 

Although (: x) :) is finite in all well behaved quantum states, S^ap^x, x') is 
singular in the limit that x and x' are null separated. There is a useful decomposition 
of operator products which classifies the various contributions to the correlation func- 
tion. Let T[x) =:<f)i{x)<f)2{x) ', where 0i and 02 are free quantum fields, or derivatives 
of fields. The stress tensor operator : : is a sum of such terms. It may be shown 
using Wick's theorem that 

T(x)T(x') = S + S 1 + S 2 , (4) 

where 

S Q = ((f) 1 (x)(f) 1 (x')) ((f)2(x)(f)2(x')) + ((pl(x)(f)2(x')) ((j)2(x)(j) 1 (x')) , (5) 



+ :<t> 2 {x)4> 1 {x'): (Mx)<hW))o+ ■M^M^)- (M^M^'))o, (6) 

and 

s 2 =:Mx)Mx)Mx')Mx')-- • (7) 

Thus the operator product T(x) T(x') consists of a purely vacuum part S , a fully 
normal-ordered part S 2 , and a part Si which is a cross term between the vacuum and 
normal-ordered parts. 
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The expectation value of the product of stress tensors can be expressed as the sum 
of these three contributions: 



{:T IJ , v {x)T a p{x'):) = {:T^ v (x) : :T a p{x') :) + (T^(x)T a(3 (x')) CTOSS 

vac i v u / 

where the last two terms are the cross term and vacuum term, respectively. In the 
limit that x' — > x, 

(T^^T^x')} cross ~ , (9) 



(x - x'^ 4 



and 



{T^(x)T a p(x')) vac ~ - . (10) 

\ t\j Jb } 



Observable quantities will be expressed as integrals of the correlation function. How- 
ever, the singularities in the cross and vacuum terms are not integrable in a naive 
way. As discussed in Refs. jHHE^ one solution to this problem is an integration by 
parts procedure. This method allows one to define the relevant integrals, and will be 
discussed in more detail below. An alternative approach, which will not be pursued 
in detail in this paper, is to use dimensional regularization. In this approach, the 
divergent part of a stress tensor correlation function can be expressed in a local form, 
which arises from counterterms in the action of the form of R 2 and R^ U R^ V , where 
R^y is the Ricci tensor [lq . 



III. RICCI TENSOR FLUCTUATIONS AND RAYCHAUDHURI'S EQUA- 
TION 



A. Raychaudhuri's Equation as a Langevin Equation 

For a classical spacetime, Einstein's equations relate the Ricci tensor R^ to the 
stress tensor by 

^ = 8tt(t^-^t;) . (11) 

In the case of a spacetime with fluctuations driven by stress tensor fluctuations, we 
assume that this algebraic relation between the Ricci and stress tensors still holds. 
A useful probe of the fluctuations of the Ricci tensor is the Raychaudhuri equation. 

Consider a congruence of either timelike or null geodesies with affine parameter A 
and tangent vector field k^. The Raychaudhuri equation gives the rate of change of 
the expansion 8 along the congruence to be |l9| 

d0 

— = -R, u k»k» -a9 2 - o^ v + . (12) 

Here a^ v is the shear and uo^ u is the vorticity of the congruence. The constant 
a = 1/2 for null geodesies, and a = 1/3 for timelike geodesies. We are interested in 
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interpreting this equation as a Langevin equation in which the Ricci tensor fluctuates. 
For the purposes of this paper, we will further assume that the shear and vorticity of 
the congruence vanishes, and that the expansion remains sufficiently small that the 
9 2 term can also be ignored, so we can write 

df) 

— = -R^k", (13) 

where is determined by the fluctuating stress tensor by Eq. (fTTj) . 
Let the Ricci tensor correlation function be defined to be 

Cpv a p(x,x') = (R fMl/ (x)R af3 (x')) - (R^(x))(R a p(x')) . (14) 

The variance of the expansion can be expressed as a double integral of this correlation 
function as 

(6 2 ) - (6) 2 = ((A6) 2 } = ( X °d\ f X °dX f C^(X,X')k^X)k"(X)k a (X')k^X / ). (15) 

Jo Jo 

Here we are assuming that the fluctuations of the expansion occur in a finite affine 
parameter interval, between A = and A = Ao- In the case of a flat background, we 
can take the stress tensor to be normal ordered. Then the Ricci tensor correlation 
function C^a^x^x') is related to the correlation function for the normal ordered 
stress tensor by 

C^a^X, X') = 647T 2 (S^ua/l ~ ^Qp, S P pa (3 ~ + -Q^ g af3 S P p ^ , (16) 

and can also be decomposed into fully normal ordered, cross, and vacuum terms. 

In the remainder of this paper, we will take the matter field to be a massless, 
minimally coupled scalar field, for which the stress tensor is 

Tp, = 0, M 0,„ - -g^ (p' p (p,p , (17) 

and the Ricci tensor is 

iV = 87r< M,- ( 18 ) 
This will now be interpreted as a normal ordered operator, so the correlation function 
Cn Ua f3{x,x') can be expressed as a sum of a fully normal ordered term, a cross term, 
and a vacuum term: 

C^(x, x') = C^(x, x') NO + C^(x, x') cross + C^{x, x') vac . (19) 

Each of these terms can be found explicitly from Eqs. ©, ©, and (jJJ) to be 

C pva \x, x') NO = 64tt 2 (: dy{x) d u (j){x) d' a <j){x') <9 ,/3 0(x') :> , (20) 

C^ af3 {x,x') cross = 64tt 2 (d p d' a D d v d'^D + d v d' a D d^d ll3 D 

+ d v & p D d p d' a D + d^d'^D d v d' a D ) , (21) 



and 

C^( X} x ') vac = 64vr 2 (d»& a D d u & p D Q + d u d ,a D d»d ,(} D ) . (22) 

Here <9 M denotes differentiation at point x and d' a that at x'. The vacuum two-point 
function is 

D = (0(x)0(x')) o = An2 x , )2 , (23) 

and 

D = (:^x)4>(x'):) (24) 

is the renormalized two-point function in the chosen state. 

The central quantity which we calculate is the variance of the expansion parameter, 
((A6) 2 ). However, it will be useful to relate it to more directly observable quanti- 
ties. Two possibilities are characteristic angle by which an images will be blurred by 
metric fluctuations, Atp, and the fractional fluctuations in the apparent luminosity 
of the source, AL/L. In Appendix A, we give a heuristic treatment of the relation 
between Acp and A9, and argue that we should expect a blurring angle of the order 
of Aip = sA6/2, where s is the distance to the source. Our treatment is heuristic 
because we actually analyze the situation of a source in flat spacetime whose posi- 
tion fluctuates. We conjecture that a careful treatment of the effects of a fluctuating 
spacetime geometry will give a result of the same order of magnitude. Fluctuations 
of 9 are also related to fluctuations in the brightness, or apparent luminosity of the 
source, as is discussed in Appendix B. In general, AL/L is given in terms of an inte- 
gral of the expansion correlation function, however, in some cases, it is of the same 
order as Aip. 



B. Averaging over a Bundle of Geodesies 

In order for the contributions of the cross and vacuum terms to be finite, it is 
necessary to average over a spacetime volume. This volume is the interior of the 
world tube defined by the bundle of geodesies. Thus in Eq. (JT5j) , we replace the 
integrations on the affine parameter by four-dimensional spacetime integrations and 
write 

((A9) 2 ) = Jd 4 xJ d 4 x" f(x) f(x') C^ va(3 (x, x') ¥ l {x)k v {x) k a {x')k p {x') . (25) 

where f(x) is the function which describes the shape of the world tube. 

It is convenient to express this averaging in null coordinates. Let the unperturbed 
rays propagate in the ^-direction, and let u = t — x and v = t + x be the null 
coordinates. We can take the affine parameter to be the v- coordinate, and then 
introduce additional averaging in the u, y, and z directions. We will take the averaging 
functions to be Lorenztian or Gaussian functions in each of these directions. Let 

g L (u,a) = —^ — (26) 
tt(u z + or) 
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and 

9G (u,a) = -^e- u2 / a2 (27) 



jit a 

so that 



/OO rOO 
gi{u,a)du = / g G (u,a)du=l. (28) 
-oo J —oo 

Equation (|23j) involves an averaging in six variables, u, u' , y, y', z, and z', as well 
as integrations over v and v'. However, in the cases that will be of interest, the 
correlation function depends only upon the differences Am = u — u' , etc. We can 
exploit the property of Lorenztian and Gaussian functions that averaging in u and u' 
separately is equivalent to an average on Am. In particular, 

/oo roo roo 

du / du' g L (u, a) g L (u', a) F(u - v!) = / dAu g L (Au, 2a) F(Au) , (29) 

-oo J—oo J —oo 

and 

/oo roo roo 

du / du' g G {u, a) g G (u', a) F(u - u') = / dAu g G {Au, V2a) F(Au) . (30) 

-oo J—oo J—oo 

Then we need only average over the differences Am, Ay, and Az. It will be convenient 
to let the averaging function be a Lorenztian of width a in Au and Gaussians of width 
b in each of Ay and Az, so 

f{x) = g L {Au, a) g G {Ay, b) g G (Az, b) . (31) 

The integrations on v and v' already amount to averaging, so we do not need to 
introduce any additional smearing in these variables. If the integrand is a function 
of the difference Av = v — v' then we may use the relation 



rvo fVQ rvo fU 

/ dv / dv 1 F(v -v')= / dAv (v - Av) F(Av) + / dAv (v + Av) F(Av) . 

JO JO Jo J-v 

(32) 

If F(Av ) is an even function, this relation becomes 

rvo rvo rvo 

/ dv / dv' F(v -v') = 2 / dAv (v - Av) F{Av) . (33) 
Jo Jo Jo 

Henceforth, we will slightly change the notation, dropping the "A"s, and write 
v, u, y, and z for Av, Au, Ay, and Az, respectively. Now our expression for the 
averaged variance in the expansion becomes 



((A6Y) 



2 



dv (v — v) + dv (v + v) 



dugL(u,a) x 



dy / dz g G (y, b) g G (z, b) C^x, x') k»(x)k"(x) k a (x')kP(x') . (34) 
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C. Analysis of the Cross and Vacuum Terms 



Here we wish to look at the state-independent and potentially singular factors in 
the cross and the vacuum terms. Let 

a = i(x-x') 2 = ^Ax" Ax^. (35) 

Then the vacuum two-point function is 

D ^ = sk- < 36 > 

Further let 

B = k,k a d»d' a D = -{W)D' - {k^fD^ . (37) 

Here D' and Dq are derivatives of D with respect to a. The contribution of the 
vacuum term to C^ va p k^k v k"^, the integrand in Eq. (jMj), is 

c^fk^KKk^ns^B 2 . (38) 

Similarly, the contribution of the cross term is 

CToSf k,k u k a k p = 256tt 2 k^k a (d»d' a D) B . (39) 

The contributions of these terms must be defined by averaging, as in Eq. (|25|). 
However, we can take the tangent vector k^ to be null, even under the averaging, so 
all of the rays in the bundle over which we average have the same direction. Then 
k^kfj, = and the first term in B vanishes. We can take k^ = —5™. then 

2m 2 

B = ~ 2T2 ' ( 40 ) 

where r 2 = y 2 + z 2 . The vacuum contribution to ((A8) 2 ) becomes 

., . n , 2 , 1024 r° , / 
((A6) 2 ) vac = — — / dv (v - v) x 

7T Z JO 



du dy dzg L (u,a) g G (y,b) g G (z,b) 



u 4 



7r 2 (r 2 — uv) 6 



(41) 



Note that once the integration on u is performed, the integrand is an even function 
of v, so we can use Eq. (J33J). 

Let us first consider the w-integration. The integrand has simple poles at u = ±ia, 
coming from the sampling function gi(u,a). However, it also has a sixth-order pole 
at u = r 2 /v, which reflects the singular nature of the vacuum term. As is well-known, 
contour integrals around higher order poles may be evaluated by integration by parts. 
Here the answer must be a real quantity, so we can write the integral on u as the 
average of two integrals, one over a contour closed in the upper-half plane and going 
above the sixth-order pole on the real axis, and a second over a contour closed in the 
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lower-half plane and going below this pole. The result is finite and just the real part 
of the residue of the pole at either u = ia or u = —ia, 



dugi(u, a) 



u 



(r 2 — uv) G 



Re 



(r 2 — iav) 6 



(42) 



It is convenient to perform the u-integration next, which yields 



rv 

/ dv (v — v) Re 
Jo 



r z — lav) 



a 4 vZ(a 6 v% + Aa\^r A + ha'vir* + lOr 



2„,2„,8 



,12\ 



20r 8 (a 2 f 2 + r 



„4^4 



(43) 



In the limit that vq ^> r 2 /a, this approaches a 2 / (20r 8 ). Both the exact form, Eq. ()43)1 . 
and the large vq limit have the same r behavior for small r. In either further 
averaging in the transverse directions is needed to produce a finite result. 
In the large limit, we can write 



256 
5^ 



a { —7 



(44) 



The evaluation of quantities such as (1/r 8 ) is discussed in Appendix C, where it is 
argued that 

where Ag has at most a logarithmic dependence on b. This leads to our final result 
for the contribution of the vacuum term in the large Vq limit, 



((A9) 2 ) vac = t P 



256 A 8 a 2 
5tt 2 6 8 



(46) 



where ip is the Planck length. 

Note that Vq = 2 s, where s is the spatial separation between the source and 
detector, and that the large Vq limit now means vq ^> b 2 /a. The fact that ((A6>) 2 ) mc 
approaches a finite value in this limit indicates that the focusing fluctuations in the 
Minkowski vacuum do not accumulate, but they can be non-zero. This raises the 
question as to whether they lead to potentially observable effects. The expected 
angular blurring or fractional luminosity fluctuations of an image at a distance s are 
given by Eqs. (jA12|) and (jB9|) . respectively, and in this case will be 



2A(p 



AL\ 



0.1 A 8 4^ 



10 



10 



-10 



cm 



10 28 cm 



(47) 



The smallest reasonable averaging scale would be of the order of the wavelength of 
the photons from the source. Even in this case, the effect predicted by Eq. (|47|) is 
too small to have been seen in any measurements that have yet been performed. 

It is also of interest to note that the asymptotic value of ((A^) 2 )^ can be made 
arbitrarily small by making the temporal averaging scale a small, although the value 
of Vq required to attain this asymptotic value grows. This behavior can be traced to 
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the factor of (fc^Ax^) 4 in Eq. (|38|). This factor in turn contributes a factor of a 4 to 
((A6>) 2 ) mc , and turns what would otherwise be an a~ 2 dependence into a 2 dependence. 
A general product of stress tensor operators needs to be averaged in time to be finite, 
but the particular combination that contribute to ((A6) 2 ) vac remains finite even in 
the limit that the time averaging scale goes to zero. 

Now let us turn to the cross term, which is determined by Eq. (|39|) . This term also 
requires averaging to render it finite. However, it is state-dependent, so we cannot 
calculate it explicitly without specifying the function D. However, we can look at 
the result of averaging the state-independent factor B, and then assuming that the 
state-dependent factor is approximately independent of v. This leads to 

((A9) 2 ) oro , B = 512n 2 k li k a (d' i ff a D) / dv(v -v)(B). (48) 

Jo 

We may evaluate the integral of (B) in the same manner as used for the vacuum 
term: 

/* o . . /*oo 2a 2 r , 9 . , r 2 (3a 2 t> 2 — r 4 ) . s 

(B}= d 2 xg(r) / dug L (u,a)B = d 2 xg(r) \ } . (49) 

J J-oo n z J (a z v z + r 4 r 



Then 



v ° , , , , „> a 2 f ,o , , a 2 fn 



. d „ K -„ )(B>= ^ j/ ^ 9(r) _^_ I _. (50) 

Note that this expression, and hence the contribution of the cross term, vanishes in 
the limit that a — > for fixed Vq. However, if we let v o — * oo for fixed, nonzero a, 
then we find 

f v 1 / 1 \ An 

dv(vo-v)(B) = -(-) = ^L, (51) 
o 7i z \r z / 7f z b z 

where Ai has at most a logarithmic dependence on b. As expected, the cross term is 
singular in the limit that b — > 0, but less so than the vacuum term. 

As was the case for the vacuum term, ((A6) 2 ) cross vanishes in the limit that a — > 0. 
In both cases, this is due to the presence of the factor of (k^Ax^) 2 in B. Note that 
if we were to take first the limit a — > and then the limit b — > 0, we would have a 
localized ray for which ((A8) 2 ) is given entirely by the fully normal ordered term. 

Here we have used a specific choice of averaging. One topic for future research 
will be to generalize this choice. In particular, it will be of interest to write covariant 
expressions for B and B 2 as derivatives of logarithmic functions, and then integrate 
over a four-dimensional spacetime region. 



D. The Physical Meaning of Smearing Scales a and b 

We have defined the singular terms of the Ricci tensor by averaging over a world 
tube containing the bundle of geodesies. The shape of this world tube is deter- 
mined by its spatial width b, which describes the transverse distance from the central 
geodesic, and its temporal width a. 
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As far as the temporal smearing is concerned, it is worth noting that the parameter 
a only appears in the pure vacuum term, and does not affect results obtained for 
state-dependent terms (the cross terms and the fully normal-ordered term). If one 
assumes that pure vacuum effects in Minkowski spacetime are unphysical, or least 
unobservable, a discussion of the meaning of a is rendered superfluous. However, 
if we remain open to the possibility that the pure vacuum effects possess physical 
significance, we should attempt to connect a with some concrete length scale in the 
problem. There are a few reasonable options to chose from. Certainly the following 
come to mind: the wavelength of the source, the coherence length of the source, and 
the duration of time over which the light source is collected. The latter, in at least 
one respect, seems the least appealing. It would suggest that the vacuum effect could 
be made arbitrarily large by leaving the measuring device on for an arbitrarily long 
time (noting the power dependence in Eq. (|47jl . If we are concerned with finding a 
upper bound when the measuring process is arbitrarily extended in time, it seems 
more likely that the coherence length of the light source, or even the smallest scale 
involved, the wavelength, should be related to the temporal smearing. The argument 
here is that the photon extends off the lightcone, but certainly no more than the 
temporal size of its wave packet, and no less than the resolution provided by its 
characteristic wavelength. 

On the other hand, it seems most natural to associate the width of the bundle 
with the parameters of the experimental set-up, determined by both conditions at 
the source and at the observer. That is, it seems we can associate the width of the 
bundle b with the two quantities: the aperture of the device used to collect the light 
rays, and the size of the emitting source. If both could be fixed to the same size, we 
would effectively enclose the bundle of geodesies within a cylinder extending along 
the entire flight path. Under such a circumstance it seems safe to say that b would 
simply be the radius of the cylinder. In the case of point source emissions, b might 
be sensibly set to the average width of the bundle along the flight path, i.e. one half 
of the aperture's diameter. For a larger source, b should perhaps an average between 
the diameter of the emitting source and the diameter of the aperture. 



IV. A THERMAL STATE 

A. Fully Normal Ordered Term in the Low Temperature Limit 

Now we turn to the case of a thermal state, so the fluctuations of the expansion 
have their origin in thermal fluctuations of the scalar field stress tensor. In this 
subsection, we consider the limit of a localized ray in the sense discussed above. The 
limit in which a and b are small is equivalent to a low temperature limit because 
the thermal wavelength is large compared to these scales. In this limit, we need to 
consider only the effect of the fully normal ordered term, so we must calculate the 
right-hand side of Eq. (|2Uj). There is a useful identity which can be used for this 
purpose: 

(: 01020304 :) = (■ 0102 :) (■ 0304 :) + (■ 0103 :) (: 0204 :) + (■ 0104 :) (: 0203 '■) , (52) 
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where the expectation values are taken in a thermal state. One proof of this identity 
is based upon the Hawking effect [2(|. First, one may use Wick's theorem to prove 
Eq. (J52j) for a generalized vacuum state ||. This is any state \ip) for which there exists 
a decomposition of the field operators into positive and negative frequency parts, 

= 0+ + 0" , (53) 

such that <f) + \ip) = ("010" = 0- However, the formation of a black hole by gravitational 
collapse transforms the in- vacuum state into a thermal state at late times |2£j] . As a 
consequence, Eq. (1H2*|) must also be valid for thermal states. 

It is also possible to give a more direct proof of this identity. Recall that the 
thermal average of an operator A can be expressed as 

w>=W- (54) 

where p is the density operator. Let us now express each of the 0j operators as a sum 
of positive and negative frequency components, as in Eq. ()5Hj) . Recall that normal 
ordering means positive frequency terms are to the right of all negative frequency 
terms. Because p has a Fock space representation as a sum of terms of the form 
|n)(n|, where \n) is a number eigenstate, only terms with equal numbers of positive 
frequency and negative frequency factors can have nonzero expectation values. Thus, 
we can write 

(: 01020304:) = (0r02 03 04 ) + (0r03 02"04 ) + (0704 02"03 ) 

+ (02 03-0101 ) + (02^04 0|03 + ) + (03 04>W) ■ (55) 

Next we can use the fact that in a thermal state 

{a\a]a k ai) = (a\ai){a]a k ) + (a\a k ) (a]ai) (56) 

to write 

(0r0 2 -0 3 f <Pi ) = (0r0l ) (02 3 + ) + (0r03 > <«4 > , (57) 

and similar relations. We can now express (: 01020304:) as a sum of products of the 
form of that on the right-hand side of Eq. ()57jl . By use of the fact that for a thermal 

(: 0102 :) = (0r0^ ) + (0 2 " 01) , (58) 

this sum can be rewritten to complete the proof of Eq. ()52|). 

By use of this identity, the fully normal ordered term, Eq. (j2*Uj) . can be expressed 

as 

C^(x, x') NO = 64vr 2 (d^d la D p d v d h D p + d u d' a D p d^d'^Dp) , (59) 

where 

D P = (: 0(*)0(*') :), = ^ M |AxP-(At + m/3)2] ' (60) 

is the renormalized thermal two-point function. Here (3 = 1/T, where T is the 
temperature, Ax = x — x', At = t — t', and the prime on the summation denotes 
that the n = term is omitted. 
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Recall that here we are discussing a limit of localized rays. Let the unperturbed 
path of the bundle of geodesies be in the +x-direction, so that Ax = At and Ay = 
Az = 0. We can take the coordinate time in a given frame of reference to be the 
affine parameter, so that A = t, X' = t' and k^ = (1, 1, 0, 0) in Eq. (fT5|). After some 
calculation, the integrand may be expressed as 



ii 



C^ vai v^v u v a v^ = C^q 1 V>" v a v 1 = -j-jp' f(k) , (61) 



where 
/(*) 



(n 2 -3k 2 )] 2 



S(n 2 + * 2 ) S J 

1 2 

= — [vr 3 fc 3 coth(vrA;) csch 2 (7r£;) + tt 2 k 2 csch 2 (irk) + irk coth(vrA;) - 3] (62) 
and k = 2 At//?. Now we can write the variance of the expansion as 

((A6) 2 ) = r dt r dt' k^ y k a k 7 = ^— r - *)/(*) . (63) 

JO JO TT Z (j° JO 

Here s is the separation of the source and detector (the flight time), k = 2s/ ft, and 
we have used Eq. (J33J). The function f(k) is finite at k = 0, /(0) = 7r 8 /8100 and 
vanishes rapidly for large k: 

/(*)~^. ( 64 ) 



<(A0) 2 > = -^, (65) 



As a consequence, both / °° dk f(k) and Jq°° dk k f(k) are finite. In this case, the large 
k limit of Eq. (|63|) is 

2 10 c 
where 

/■oo 

c = / dkfik) « 0.3468. (66) 

JO 

We can now write the root-mean-squared fluctuation in the expansion as 



A9 rms = 128 v^° p . (67) 

7T 

The fractional luminosity fluctuations of the image of a point source viewed at a 
distance s through a bath of scalar radiation at temperature T becomes, for localized 

rays, 

ff) J^e p ^. (68) 

V Li J rms V O 7T 
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Equation ({SB} can be expressed as 



(irL = °' 02 ( tpU 1 ( i^f) 3 = 10 " ( 10^) 1 (k£v) 3 ■ < 69) 

Thus a source seen at a cosmological distance, 10 28 cm « 3Gpc through a thermal 
bath at a temperature of lO 6 ^ could in principle show large luminosity fluctuations 
if one could select rays localized on a scale less than the thermal wavelength (3. Of 
course, there cannot be such a bath in our present universe, as its energy density 
would have long dominated the expansion of the universe. Note that 10 6 km is the 
Schwarzschild radius of a black hole with a mass of 3 x 1O 5 M , in the mass range of the 
supermassive black holes which may be in the centers of many galaxies, and possible 
linked to the phenomenon of gamma ray bursts. In at least some cases [22j |. these 
bursts have produced gamma rays with energies well in excess of lGeV. Although 
the details of the size and nature of the emitting regions is unclear, this raises the 
possibility that within the burst sources, there are quantum metric fluctuations large 
enough to achieve AL/L w 1 for a bundle of rays localized on a scale small compared 
to thermal wavelength (3. Any realistic observation is likely to average over a much 
larger scale, so we now turn to the opposite limit. 



B. High Temperature Limit 



In the limit in which the averaging scales a and b are large compared to the thermal 
wavelength, we need to adopt a different approach. Now the effect of the averaging 
cannot be ignored, even in the fully normal-ordered term. It is easiest to include the 
effects of the vacuum, cross, and fully normal-ordered terms all together by using the 
full thermal Green's function, 



D Fp = (:<P{x)<P{xi):) p = £ 



1 



3 4vr 2 [|Ax| 2 - (At + in(3) 2 } ' 
If we impose the same type of averaging as used in Sect. 1111 CI then we can write 
((A6f) = 



(70) 



512 



7T 



dv (v —v)+ dv (v o + v) 



I'D 



dug L (Au,a) / d 2 xg(r) (d 2 v D F/3 ) 2 .(71) 



In the high temperature (small (3) limit, we can replace the sum on n m. Dpp by an 
integral. Let £ = (3n, and let 



E 



i 



(72) 



Then we can write 



dlD 



16 



Fp 



(i£ + u) 



2r 2 — (v — u) 



tt 2 /3 J-oo * {{v - u) 2 + 4r 2 - (2i£ + v + u 



213 



2ti(3 [( w _ M )2 + 4 r 2]| 
(73) 
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This leads to 

•z 

/(A^ 2 \ = _ 

IP 



((AO) 2 ) = § / dug L (Au,a) I d 2 xg(r) x 



VQ 

dv (vq —v) + I dv (vq + v) 



[2r 2 -(v- u) 2 ] 



212 



[(v - u) 2 + 4r 2 ] 5 

In the limit that v$ is large compared to either a or b, we can write this as 

32v r°° , /A x /• , 2 , x Z" 00 , (2r 2 -v 2 ) 2 27ttw /l 



(74) 



((A0) 2 ) = — p. J^dug L (Au,a) Jd 2 xg(r) f_ 



f3 2 J-oo ' 'J " y J J-oo (f 2 + 4r 2 ) 5 512/3 2 \r 5 , 

(75) 

Note that in this limit, the integrand of the f-integration becomes independent of u, 
and hence the result is independent of a. 

The final expression involves the transverse average of a singular function, 1/r 5 . 
However, unlike the singular functions encountered earlier in the vacuum and cross 
terms, this is simply an artifact of our approximations. The dominant contribution 
here is from the fully-normal ordered term, and should be finite. The apparent 
singularity at r = arises because our approximation of replacing the sum by an 
integral assumes that f3 is small compared to all other length scales including r. An 
exact calculation of the fully-normal ordered contribution to ((AO) 2 ) would, in the 
high temperature limit, replace Eq. (|7H|) by 

(m 2 ) = ^ </(r)> , (76) 

where f(r) is a function which falls as r~ 5 for r ^ /?, but approaches a constant for 
r ;$ (3. We can estimate (f(r)) in the high temperature limit by simply cutting off 
the integration on r at a lower limit of order (3. Thus, 

(fir))* - r-e^«— (77) 

This leads to an estimate for ((AO) 2 ) of 

(m 2 ) = -^E°, (78) 

where K is a constant of order one. 

The corresponding fractional luminosity fluctuation is 

AL\ Ksl n9 

bpl P K ' 

Note that as compared to Eq. (|p^j). the effect is now suppressed by a factor of (5/b. 
This is essentially the effect of classical averaging of the luminosity fluctuations; when 
one averages over a larger source, the fluctuations are suppressed. In Eq. (J7HJ), we 
might take b to be either of the order of the size of the source, or of the average 
of the source size and the observer's aperture size. In either case, the luminosity 
fluctuations will be strongly suppressed compared to the estimates given in Eq. 
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V. DISCUSSION 



In this paper, we have treated the effects of the fluctuations of a massless, mini- 
mally coupled scalar field stress tensor on a flat background. In particular, we used 
the fluctuations in expansion of a congruence of null rays as an operational measure 
of the effects of stress tensor fluctuations on the spacetime geometry. Although the 
cross and vacuum terms are potentially singular parts of the stress tensor correlation 
function, we argued that their contributions can be rendered finite by an integration 
by parts procedure. In the limit of localized rays, they actually yield a zero contri- 
bution to the variance of the expansion. This is in contrast to other situations 0,0], 
where the cross term often yields the dominant contribution. We calculated an ex- 
plicit formulas for the fluctuations in the expansion produced by a thermal bath. This 
in turn led to a results, Eqs. (J55J) and (j79|) . for the fractional luminosity fluctuations 
of the image of an object seen through a fluctuating spacetime filled with a thermal 
bath. Although these effects are usually small, there are conceivable astrophysical 
situations where AL/L would become large if the rays could be sufficiently localized. 

The effects of the fully-normal ordered term are unambiguous, but those of the 
vacuum and cross terms contain singular parts which need to be defined by a space- 
time averaging over a bundle of rays. This may be done by the integration by parts 
procedure discussed in Appendix C. This procedure involves some undetermined con- 
stants, some of which can be eliminated by a scaling argument. Unfortunately, the 
overall constant which determines the magnitude of the vacuum contribution is not 
fixed. Further work is needed to determine this constant uniquely. As noted above, 
one topic for future research will be the use of covariant averaging. 

Note that Eq. (ffiTj) describes the effects of the Ricci tensor fluctuations, not the 
expansion due to any mean Ricci curvature. Our treatment has assumed a nearly flat 
background spacetime, More generally, the non-zero mean value of the stress tensor 
will lead to mean Ricci curvature and hence classical focusing. 

Recall that we assumed that the 2 term in the Raychaudhuri equation, Eq. ((12} , 
can be neglected. As ((A6) 2 ) grows, this assumption will eventually break down. 
However, in many cases ((A6) 2 ) remains small compared to the Ricci tensor term 
for a very long time. We can give an explicit estimate of this time for the case of 
the thermal bath. In this case, the root-mean-square fluctuations in the Ricci tensor 
term are of order {R^v^v v \ rms ~ fpT 4 . Thus {{AO) 2 ) <C \R^ LV v ll v u \ rms so long as 

s « ^ ■ (80) 

For temperatures less than about 1 GeV, this condition will be satisfied for any s 
less than the size of the observable universe. We also assumed that the shear and 
vorticity terms in the Raychaudhuri equation can be ignored. So long as the bundle 
of geodesies is hypersurface orthogonal, the vorticity remains zero in all classical 
spacetimes, so it is reasonable to expect that small quantum geometry fluctuations 
will not generate any vorticity. The shear can grow in response to non-zero Weyl 
curvature. Thus if the stress tensor fluctuations cause the Weyl tensor to fluctuate, 
through the excitation of gravitational wave modes, then eventually {a^ u a^ v ) could 
become important. However, it seems reasonable to speculate that it will not grow 
any more rapidly than does {{AO) 2 ). 
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In summary, we have illustrated how expansion fluctuations provide a concrete 
measure of the effects of passive quantum geometry fluctuations. 
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APPENDIX A 

In this Appendix, we will illustrate the geometric meaning of the expansion 9 
in an explicit example {HI], and then use this example to infer a relation between 
fluctuations of the expansion and the angular blurring of an image. First consider a 
point source in flat spacetime which is emitting null rays, whose tangent vector field is 
f M = dx^/dX. Choose a coordinate system in which the origin of spatial coordinates 
is at the location of the point source, and take the affine parameter to be coordinate 
time in this frame: A = t. Then 

w" = (l,v), (Al) 

where v is a unit 3-vector, which is also the unit position vector in this coordinate 
system 

v = x = ^ . (A2) 
|x| 

The expansion is defined as the trace of the tensor v ^ in the subspace orthogonal 
to v**: 

6 = Vv»(g ia , + v ll v v ). (A3) 
In our case, the tensor V u has only spatial components 



V J v l = & v l = — ^ . (A4) 

l x l 

Now 9 can be found to be 

9 = &> v l (<% + viOj) = ~ , (A5) 

where s = |x| is the distance from the source. If we repeat this calculation for the 

case of timelike geodesies, using the 4-velocity w M as the tangent vector, the result is 

«=^, (A6) 
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FIG. 1: A source S at a mean distance s from an observer O is displaced in direction tp by 
a distance p from its average position A. The effect is to shift the apparent position of the 
source by an angle A(p. 



where 7 = 1/y/l — v 2 . In both cases, 9 oc 1/s, so as the distance from the source 
increases and the rays become more nearly parallel, the expansion decreases. 

Now let us relate the fluctuations of 9 to fluctuations in the apparent angular 
position of the source. In our case, the ^-fluctuations are due to fluctuations in 
the spacetime geometry between the source and the observer. However, one could 
achieve the same effect with a source whose position fluctuates in a fixed Minkowski 
spacetime. Consider the situation illustrated in Fig. Q where the source is displaced 
from its mean position by a distance p in a direction given by the angle if). For an 
observer at a mean distance of s, the expansion changes by 

2 2 

A9 = — - As = — p cosip, (A7) 
and the angular position of the source changes by approximately 

A V =^, (A8) 
s 

provided that p <^ s. We now assume that p and if) fluctuate independently of one 
another and that all values of if) are equally probable. Then 

(p 2 sin 2 ^) = (p 2 )(sin 2 ^) = i(p 2 ) (A9) 

and 

(p 2 cos 2 ^) = (p 2 )(cos 2 ^) = i(p 2 ). (AlO) 



18 



This would happen, for example, if p does not fluctuate and ip fluctuates randomly. 
This allow us to relate the expansion fluctuations to the angular position fluctuations: 

{(A^) = \s 2 ((Afff) . (All) 

Thus the effect of ^-fluctuations is to blur the angular resolution of the image of the 
source by an amount 

Ay? rms = X - s A6 rms . (A12) 

Strictly, we have treated only a source in flat spacetime whose position fluctuates. 
However, we offer this as a heuristic estimate of the magnitude of angular blurring 
expected in a fluctuating spacetime. 



APPENDIX B 

In this appendix, we show how the stochastic behavior of the expansion leads to 
small deviations in the density of geodesies. In the case of light rays this translates 
into luminosity fluctuations. The connection between luminosity fluctuations and 
expansion fluctuations can be seen as follows. Let us write the expansion as the sum 
of a classical (deterministic) component 9 C and a stochastic component 9: 

6(t) = 6 c (t)+6(t). (Bl) 

Here t is an affine parameter along the bundle of geodesies. The stochastic component 
can expresse as 

at 

where A is the cross sectional area of a bundle. Now we assume 9 is a small stochastic 
perturbation. This means 

A = A + A c w Ac, (B3) 

where A is the stochastic variation in the area (zero at t — 0). Using the Taylor 
expansion of log(l+x) for small x, we find 



/ 



d(A/A c ) 

e « -L_J ( B4) 

* A Ao 

dt' 6(t') « 4- = — , (B5) 

' : A C P 



where p is the density of geodesies. Since the number of geodesies, n, is conserved 
(total energy flow), while the transverse area A is distorted (density fluctuations), 

Ti 

Intensity or Luminosity = —. (B6) 
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The fractional luminosity fluctuation is therefore given by 

AL Ap rt 



« / at' eit'). (B7) 

Jtn 



L p Jto 

Now the variance of the fractional luminosity fluctuations can be computed as 

( (^~) 2 ) = f r dt> dt " { m e{n } ■ m 

For large flight times (where s in the upper bound of the integral is much greater 
than the correlation length scale), there is a simple connection between the variance 
of the fractional luminosity fluctuations and the variance of the expansion: 

((^) 2 )- * 2 (^)}- (B9) 

Let us show this for a typical positive-definite correlation function, Ci(At). We 
have 

( 6(t) 0(f) ) = fdtx f dt 2 Ci(|At|). (BIO) 
Jo Jo 

The function C\ is symmetric in its arguments, so just consider t < t', then sym- 
metrize later. The first integration over dt 2 produces the area under the Ricci cor- 
relation curve, c . (We will be considering large flight times, so the edge effect at 
the origin is neglected.) After the second integration, and symmetrizing, we have the 
expression for the expansion correlation function: 

( S{t) 6(t') ) = c t, t' > t, (Bll) 
= c t', t > t'. (B12) 

Now we can calculate the fractional luminosity fluctuations: 



f f dt dt' ( 9{t) 9(f) ) 
Jo Jo 


(B13) 


2 f dt f dt' ( 6{t) 0(t') ) , t>t', 
Jo Jo 

2 c f S dt ( dt' t' 
Jo Jo 


(B14) 
(B15) 


c / dt t 2 
Jo 


(B16) 


1 3 


(B17) 


^ 2 (^) 2 )- 


(B18) 



Explicit calculations show that Eq. (jB9|) also holds for the vacuum correlation 
function encountered in this paper. In the case of the vacuum term, direct calculation 
reveals 

((f) 2 ) =i«"<w"L- < B19 > 
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APPENDIX C 



In this Appendix, we will discuss two methods for defining quantities such as 
(1/r 8 ). The first approach involves an integration by parts. Let us first illustrate this 
approach in one space dimension with the quantity 



P 



(i> 



dxx n gi(x) 



(CI) 



where gi(x) = e x2 ^/(y/ir b). If n > —1, then the above integral is absolutely 
convergent and there is no ambiguity in the definition of . If, however, n < — 1, 
then the integral in Eq. (jClj) diverges at x — 0. In the case that n is a negative 
integer, let m = —n and write 



\m— 1 



d r > 



X 



X" 



2 (m - 1)! dx r 



In 



x 



■Pm-l(^) 



(C2) 



where Xq is a constant and P m _i(x) is an arbitrary (m — l)-th degree polynomial in 
x. We can now insert this identity into Eq. (|C1|) and integrate by parts to define 



rn 



X' 



2(m- 1)! 



dx 



In 



x 



+ P m -i(x) 



d r ' 



dx r 



■9i(x). (C3) 



The resulting integrals are not only convergent, but are also independent of the 
choices of xq and the coefficients of P TO _i(x). (These are not really independent, as a 
shift in xq shifts the zeroth order term in P m _i.) Furthermore, the results are always 
of order b~ m . For example, 

(l/x 2 } = -2/b\ (C4) 
(1/x 4 ) = 4/(36 4 ) , (C5) 
(1/a; 6 ) = -8/(156 6 ) , (C6) 

and 

(1/x 8 ) = 16/(1056 8 ) . (C7) 

Note that if we were to allow the polynomial P m -i to be a function of \x\, rather 
than x, then the final result will depend upon the coefficients of the odd powers of 
\x\. For example, let 



Pi = Pi(\x\) = c + ci |x| . 



Then we find 



P-2 



2 

&2 



ix b 



(C8) 



(C9) 



In effect, what has happened here is the addition of a term proportional to S(x) into 

the representation of 1/x 2 , which has caused p_ 2 to acquire a term proportional to 
<?i(0) with an arbitrary coefficient. 
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The second approach to defining the p^ involves a recurrence relation. Note that 
if we let a = 1/b 2 , then we can write 



dgi 
da 



This allows us to derive the recurrence relation 



o {1) 

Pn+2 



da 



a 



(CIO) 



(Cll) 



We can start from p^ = 1 and uniquely generate all p^\ where n is a positive even 
integer from Eq. (|C11|) . We can also integrate this relation to define p^ for negative 
even values of n, but the result will contain undetermined constants of integration 
which are equivalent to the arbitrary coefficients in P m _i(|x|). 

Now we turn to two space dimensions and examine the quantities 



P (2) 

rn 



d x r n g{r) = 2ir / dr r n+ g (r ) 



(C12) 



where the sampling function is now g(r) = e~ r ' /(irb ). Here the integrals converge 
if n > —2, but are divergent for n < —2. In the latter case, we can still use integration 
by parts to define pffl. Consider first the case where n = —2£ is a negative, even 
integer, and write 



1 



(V 2 



where 



In 2 



1 d 



+ ft(r) 



V - - I ri- 



r dr \ dr / 

is the Laplacian operator and f m {f) is any function which satisfies 



(VY ft(r) 







(C13) 
(C14) 

(C15) 



for r ^ 0. If we insert Eq. (jC13|) into Eq. (jC12|) . the Laplacian operators may be 
successively moved to act on g(r) by repeated use of Green's theorem: 



J {FV 2 G - GV 2 F)dV = j(FVG - GVF)dc 



(C16) 



Here the region of integration is a disk of radius R. In the limit that R — > oo, the 
surface term vanishes because g(r) and its derivatives vanish rapidly. The result is 



o {2) 

P-21 



dr r 



a e In 2 U> + fi(r) 



iy 2 Yg(r) 



(C17) 



This integral is convergent, but depends upon various arbitrary constants. 



22 



This is best illustrated with some specific examples. First let £ = 1, in which case 
we find ai = 1/8 and 

' r 



fi = c In 



ci 



This leads to the result 



1 

&2 



2 In 



+ 7 + c 



(Cll 



(C19) 



where 7 is Euler's constant. There is really one arbitrary constant in this result, as 
Co may be absorbed into a redefinition of the constant r . We can understand the 
appearance of c in Eq. (jC19|) from the fact that 



VVi(r-) = -27rco<5(x) 



(C20) 



in two dimensions. Similar results arise for larger values of £, but with more constants, 
For £ = 2, where a<i = 1/32, we can solve V 2 /2 = fi to find 



/2 = iCo r 2 



In 



- 1 



+ c 2 In 



and 



1 



1 

1? 



21n(y) + 7 + l-8c 2 



r 
+ 



For £ = 3, a 3 = 1/512 and 
1 



/a 



128 
which yields 



c r 



2 In ( - 

Cl 



1 

456 



r 



In ( - 

C3 



4 In ( y ) + 2 7 + 3 - 64c 4 



Finally for £ = 4, where a 4 = 1/18432, 
1 



13824 

1 2 

4 C4r 



c r 



6 In 



/ r 

VCi 



In ( - 

c 3 



11 



c 6 In — 

^c 7 



H Co r 

128 



2cq 
b 2 



c 4 In ( — 



8c 2 2c 

b 4 + & 2 



2 in ( - 

Cl 



(C21) 
(C22) 

(C23) 
(C24) 



(C25) 



and 



366 8 



12 In ( j ) + (J- I I 



768c 6 64c 4 _ 8c 2 2c 



IT + TF • ( C26 ) 



Now let us consider the case where n = —2£—l is a negative odd integer and write 



1 



(V 2 



\(2£-l] 



H12, 



+ ft(r) 



(C27) 
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A procedure analogous to that in the case of even n leads to 



[(2€-l)!!] 



+ rf e (r) 



(V 2 ) £ ^(r). 



The results for 



1 and 



2 are 



r 3 / 6 3 6 2 



and 



1 



4^ 



8c 2 2c 
6 4 6 2 



We can also derive a recurrence relation for p^ 

dg_ = n 

da \a 



Note that 



9- 



This leads to the relation 



n {2) 

Pn+2 



da \ a J 



(C28) 



(C29) 



(C30) 



(C31) 



(C32) 



We may start from either p^ 1 = 1 or p^[ = ^fn/b and integrate this recurrence 
relation to find pffi for all negative integers n. The results agree with those of the 
integration by parts approach, with the constants q arising as constants of integra- 
tion. 

Now we wish to address the difficult question of whether these constants can be 
determined by physical arguments. Let us first consider the one-dimensional case. 
So long as we require the polynomial P m -\ to be a function of x, as opposed to \x\, 
then 

d m 

— P m _ 1 = 0, (C33) 

and there is no ambiguity. The undetermined constants which appear if P m ~i is a 
polynomial in \x\, or if we integrate the recurrence relation, all have dimensions of 
inverse powers of length. One way to resolve the ambiguity is to require that p^ 
be proportional to b n , which sets all of the constants of integration to zero. This is 
essentially a scaling requirement, that if b — * 2 b, for example, then pff — > 2 n p^\ 
Another way to say this is that p$ should depend only upon the sampling scale b 
and not on any additional length scales. 

For odd n, we can apply the same requirement to the two-dimensional case, and 
similarly require that 

pi 2) oc b n . (C34) 

In the case that n is a negative odd integer, this uniquely defines p^ by setting the 
constants C\ = 0. For example, now 



2y^F 

b 3 



(C35) 
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and 




(C36) 



For n even and negative, there is still an undetermined multiplicative constant asso- 
ciated with a term logarithmic in b, so we will have 

<^) = w- (C37) 

where we allow A 2 e to possess a possible logarithmic dependence on b. 
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